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ABSTRACT 

We have studied the variability of PSR J0529— 6652, a radio pulsar in the 
LMC, using observations conducted at 1390 MHz with the Parkes 64-m tele- 
scope. PSR J0529— 6652 is detectable as a single pulse emitter, with amplitudes 
that classify the pulses as giant pulses. This makes PSR J0529— 6652 the second 
known giant pulse emitter in the LMC, after PSR B0540— 69. The fraction of 
the emitted pulses detectable from PSR J0529— 6652 at this frequency is roughly 
two orders of magnitude greater than it is for either PSR B0540— 69 or the Crab 
pulsar (if the latter were located in the LMC). We have measured a pulse nulling 
fraction of 83.3 ± 1.5% and an intrinsic modulation index of 4.07 ± 0.29 for PSR 
J0529— 6652. The modulation index is significantly larger than values previ- 
ously measured for typical radio pulsars but is comparable to values reported 
for members of several other neutron star classes. The large modulation index, 
giant pulses, and large nulling fraction suggest that this pulsar is phenomeno- 
logically more similar to these other, more variable sources, despite having spin 
and physical characteristics that are typical of the unrecycled radio pulsar pop- 
ulation. The large modulation index also does not appear to be consistent with 
the smal l value predict e d for this pulsar by a model of polar cap emission out- 



lined by iGil fc Sendykl (120001 ) . This conclusion depends to some extent on the 
assumption that PSR J0529— 6652 is exhibiting core emission, as suggested by its 
simple profile morphology, narrow profile width, and previously measured profile 
polarization characteristics. 

Subject headings: pulsars: individual (PSR J0529— 6652) - radiation mechanisms: 
non-thermal 
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1. Introduction 

From the time of their discovery more than 40 years ago, the detailed mechanism by 
which pulsars generate radio emission has remained elusive. Models have been proposed over 



that time to explain pulsar emission, and in principle these mode 
observations. One relatively recent emissio n model, proposed by 



s can be tested with pulsar 



Gil fc Sendvkl feOOOf ). is an 



extension of the polar cap spark model of Ruderman fc Sutherland ( 19751 ) in which bursts 
of plasma ("sparks") are produced at the polar cap. In this model, instabilities in the 
plasma traveling along the magnetic field lines induce radio emission. The model postulates 
a densely packed region of these sparks, with a characteristic size and separation h. One 
can determine a characteristic number of sparks a (in one dimension) across the polar cap 
according to 



r p /h, 



(1) 



where r p is the polar cap radius. In the iGil fc Sendvkl (120001 ) model, a is a "complexity 
parameter" that can also be determined observationally by 



5(P/10- 15 ) 2/7 (P/1 s) 
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(2) 



where P and P are the pulsar's spin period and period derivative, respectively. There are 
several ot her similar models having d ifferent scalings to determine the complexity parameter 
(see, e.g.. iBurke-Spolaor et al. ( 2012 ) for a brief summary and references to these models). 
Jenet &: Gill (120031 ) provide details on how such models might be tested by comparing a 
pulsar's observed intensity modulat ion with the variabili ty predicted by the complexity pa- 
rameter in the model. Although the lGil fc Sendyk (120001) mode l does not make a distinction 
between core and conal emission from the pulsar (jRankinlll983l ). a model test is most useful 
using pulsars with core emission, since observationally it is the most straightforward probe 
of the polar cap region, it is not affected by viewing geometry effects from an off-center cut 
through the emission region, and it is less likel y to exhibit cone-based p ulse drifting effects 
which would increase the observed modulation (IWeltevrede et al.ll2006al ). 



One measure of a pulsar's intensity variability that can be employed in such tests is the 
phase-resolved modulation index, m(0). Here is the longitude of the pulse phase, and m(0) 
is the ratio of the standard deviation <r(0) of the observ ed pulse intensities to the mean pulse 



intens ity averaged over many adjacent pulses, (/(0)) (IJenet fc Gilll2003l ; IWeltevrede et al. 
2006ah : 
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As outlined in lJenet fc Gill (120031 ). one would expect an anti-correlation between a and 
m owing to the likelihood that overlapping sparks in a complex emission region (large a) 
would wash out and tend to reduce the degree of modulation, creating more steady emission 
(smaller m). The phase information is important to retain since the measured modulation 
index near the edges of a profile peak usually increases relative to the center, in some cases 
because conal emission may be present. The center of the peak of the pulse profile therefore 
offers the best measurement of the properties of the core emission. The largest signal is also 
present at the profile peak, so the most precise constraint on m(4>) can generally be made 
at that pulse phase. 

The modulat i on ind ex has previously been measured for a number of radio pulsars. 
Weltevrede et al.l ( l2006al ) observed 187 pulsars at 1400 MHz to investigate the modula- 
tion and sub-pulse drifting properties of the sa mple, and they produced m odulation index 
measurements for 175 of them. More recently, iBurke-Spolaor et al.l (120121) reported mea - 
surements fo r a sam ple of 103 pulsars from the HTRU 1400-MHz survey (IKeith et al.ll2010l ). 
Jenet &: Gill (120031 ) analyzed 2 pulsars observ ed near 1400 MHz and re-analyzed an other 
10 that were previ ously observed at 430 MHz ( IJenet et al.ll2001t IWeisberg et al.lll986l ) . The 
Jenet fc Gill (120031) sample includ es a subset of 28 pulsars analyzed by IWeisberg et al.l (119861 ) 
at 430 MHz. I Jenet fc Gill ( 120041 ) also measured m for PSR B1937+21 at 1400 MHz. Other 



measu rements of m (but not phase resolved) were reported in Table 4 of IWeltevrede et al. 



( 1201 it ) for members of several neutron star classes having differ ent types of emission. These 
objects were the Crab pulsar - a classical giant pulse emitter ( Staelin fc Reifenstein Il968 



Argyle fc Gowerlll972l ). the radio-emitting magnet ar XTE J1810— 197 dCamilo et al. 



ble RRAT-link pulsar PSR B0656+14 ( IWeltevrede et al. 



2006|), 



the rotating radio transient (RRAT) J1 819— 1458 (McLaughlin et al. 20061 ). and the possi- 



2006bl ). The modulation indices 



reported for three of these four sources were for frequencies near 1400 MHz; PSR B0656+14 
was measured at 327 MHz. 

Generally, m is observed to be small, with typica l values between about 0.5 and 1 
( IWeltevrede et al.ll2006al 1201 It IBurke-Spolaor et al.ll2012l). T he histogram of measured mod- 



Weltevrede et al. J2006a|)_ shows no measured m with a 



ulation indices shown in Fig. 8 of 
value above 2, and only two of the 103 measure ments of IBurke-Spolaor et al.l (120121) hay e 
m > 2. As seen in Fig. 1 of I Jenet fc Gil fj2003h and Fig. 10 of IWeltevrede et al.l (bo06ah . 
the anti-co rrelation between m and a i n the observed sample of pulsars is weak. This is also 
seen in the IBurke-Spolaor et al.l ( 120121 ) sample. One reason may be the inclusion of pulsars 
with conal emission in these larger samples; the modulation index and complexity parameter 
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are not predic ted to be as strongly correlated for conal emi tters as they are for pulsars with 
core emission (IJenet &: Gilll2003l ; iBurke-Spolaor et al.ll2012l ). Another difficulty in establish- 
ing correlations is that the range of m measured in the observed data set is small. One 
well-known exception to this is the Crab pulsar. The Crab's modulation index is due to its 



well-studied giant pulses (iStae 



in Sz Reifenstein 



1968 



Argyle fc Gowerlll972l ; iLundgren et al. 



1995c IWeltevrede et al.ll2006at iKaruppusamy et al.ll2010l ). However, this pulsar is somewhat 
atypical since it is young and very energetic and is also a well-known X-r ay emitt er. In 



addition, the Crab pulsar is difficult to classify as a core or conal emitter (iRankinl Il990l ). 
The Crab may have conal emission and other complicating geometric profile features that 
could produce a larg er measured modulation index, making its use for a simple test of the 
Gil fc Sendvkl (boOOT model difficult. 



However, PSR J052 9— 6652 is one pulsar that could be used for a consistency check of 
the iGil fc Sendvkl (120001 ) model. This is a 976-ms pulsar that was first discovered in a 600 
MHz pulsar survey of the Large Magellan ic Cloud (LMC), and it was the first extragalactic 



pulsar discovered (IMcCulloch et a l. 1983 ). Af t er its disco very, the pu l sar w as studied at 



frequencies near 600 MHz by iMcConnell et al.l (11991) and ICosta et 



results for the pulsar were subsequently published by I Crawford et al.l (120011 ) . 



(1991), and timing 



PSR J0529— 6652 was redetected as a known p ulsar in a more recent 1400 MHz Parkes 



surve y of the Magellanic Clouds (MCs) for pulsars (jCrawford et al.ll200ll ; Manchester et al. 
20061 ). In these redetections, it was evident that the pulsar was variable on a time-scale 



of minutes or less and showed possible nulling behavior, unlike any of the other MC radio 
pulsars detected. PSR J0529— 6652 also had clearly detectable single radio pulses in these 
observations, making this pulsar phenomenologically unlike every other extragalactic radio 
pulsar yet discovered^] PSR J0529— 6652 is very luminous at 1400 MHz, ranking in the top 



2% o f the known radio pulsar population in the ATNF pulsar catalog ([Manchester et al. 
20051 ). However, apart from its luminosity and apparent nulling behavior, the pulsar is 



unremarkable, with a period and period derivative that are typical of the unrecycled radio 
pulsar population. 



Previous measurements of PSR J0529— 6652 taken with Parkes at 600 MHz by I Cost a et al. 



(119911 ) showed that the pulsar has a simple, single-peaked morphology, with little pulse broad- 
ening . The width of the measured profile is ~ 3% of the period at 600 MHz (ICosta et al. 



19911 ) . The integrated pulse profile of PSR J0529-6652 at 1400 MHz is also single-peaked 



3 The LMC pulsar PSR B0540— 69 also has detectable single radio pulses (| Johnston fc Romanill2003l ) . but 
i t is not detectable as a periodic source at 1400 MHz. It is also a young, rapidly rotating X-ray emitter 
(jSeward et al.lll984l ). making it different from PSR J0529— 6652 in this respect as well. 
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and uncomplicated (see our analysis below), and its width at both frequencies is consis- 
tent with the empirical rela tion for the pulse width at 1 GHz of W = 2.45°P~ 1 / 2 / sin a 
presented by iRankinl ( 119901 ) for core emitters, if a reasonable magnetic inclination angle 
is assumed (a ~ 15°). These features are consistent with a co re emission i nterpr etation 
for PSR J0529-6652. A polarization profile at 600 MHz from ICosta et al.l jl99lh shows 
significant linear polarization (~ 25%), with weaker circular polarization that changes sign 
symmetrically near the profile peak. Thi s sign reversa l also supports the interpretation of 
PSR J0529— 6652 as a single core emitter (jRankinlll990l ). In addition, the majority of known 
pulsars (~ 70%; Rankin 1990) that have been classified have core components. On this basis, 
we conclude tha t PSR J0529— 6652 i s likely to be exhibiting core emission, making it useful 
for a test of the bil fc Sendvkl J2OO0I ) model. 



2. Observations and Analysis 



In order to investigate the pulse variability characteristics of PSR J0529— 6652, we 
observed the pulsar with the Parkes 64-m radio telescope in two separate observations of 
length 5.0 and 3.8 hr in Feb 2008. Both observations were conducted at a center frequency 
of 1390 MHz using the cent e r beam of the multibeam receiver and analog filterbank system 
( jStaveley-Smith et al.l Il996t Manchester et al.l 1200 ll ). A bandwidth of 256 MHz was split 
into 512 x 0.5 MHz channels, and each channel was 1-bit digitized and sampled every 500 
/is (see Tabled]). The observing setup was identical to the one used i n a radio pulsar searc h 
of the X-ray binary XTE J0103-728 in the Small Magellanic Cloud JCrawford et al.ll2009h . 
which was part of the same observing campaign. Radio frequency interference (RFI) can be 
a significant problem at 1400 MHz at Parkes, and much of the data for PSR J0529— 6652 
was at least partially corrupted by RFI. We selected a portion of the first observation that 
was clean of RFI for the analysis. This subset consisted of 4195.5 seconds of integration, 
corresponding to 4299 complete pulses. 

We performed the following operations on PSR J0529— 6652 as well as on three bright 
test pulsars, PSRs J0437— 4715, J0536— 7543, and J1359— 6038, in order to test our process- 
ing algorithms^ We dedispersed the raw channelized data for each pulsar using th e catalog 
dispersion measure (DM) from the ATNF pulsar catalog jManchester et al.l [200 5)FI to pro- 



4 All three test pulsars were observed with the same Parkes observing system as PSR J0529— 6652 (but 
at different epochs). In the case of PSR J1359— 6038, a slightly different center frequenc y (1374 MHz) and 
filterbank system (consisting of 96 x 3 MHz channels) were used ( Manchester et al. 200lh (see Table [lj. 

D http: / /www. atnf.csiro.au/research/pulsar/psrcat/ 
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duce a dedispersed time series. We then created a pulse stack (e.g., Weltevrede et al. 2006a; 
Burke-Spolaor et al. 2012) by folding the dedispersed time series modulo the topocentric 
pulse period and stacking the resulting consecutive pulses. This produced a two-dimensional 
array consisting of pulse number vs. pulse phase bin. In each case, 128 pulse phase bins 
were used for the pulse stack. We determined the mean pulse profile by summing the pulses 
in the pulse stack while preserving pulse phase. All elements in the pulse stack were then 
subsequently adjusted (normalized) so that the calculated mean intensity had zero mean for 
the off-pulse bins and unity peak value. The mean intensity and mean intensity squared were 
then recomputed for each bin using this normalized pulse stack. Fig. [1] shows the normalized 
pulse stack for J0529— 6652 and the three test pulsars for comparison, and Fig. [2] shows the 
mean intensity profile for PSR J0529— 6652. As is the case at 600 MHz, the mean pulse 
profile for PSR J0529— 6652 at 1390 MHz has a simple morphology and narrow width, with 
a duty cycle of 3-4 bins (corresp onding to ~ 25 ms , or ~ 3% of the pulse period). This simi- 
larity to the width measured by lCosta et al.l (1199ll ) at 600 MHz indicates that the pulse does 
not experience significant broadening at lower frequencies and remains narrow across a range 
of frequencies. The phase-resolved modulation index was measured for each profile phase 
bin using the normalized pulse stack according to Eq. El An uncertainty in each modulation 
index value was computed using the rms values of the intensity and intensity squared within 
each bin and the number of points that were summed within each bin. The modulation 
index values were corrected to account for contributions from the interstellar medium (ISM) 
(this is discussed in Section 3.2), and the resulting phase-resolved modulation index for PSR 
J0529— 6652 is plotted and discussed later in the paper. 

Nulling fractions (NFs) were also computed for PSR J0529— 6652 and the three test 
pulsars using on and off-pulse in tensity histograms that were created from the pulse stacks, 
as outlined by lWang et al.l (120071 ) . Prior to this, however, we tested ou r procedure by reco m- 
puting the NF for three known nulling pulsars measured at Parkes by IWang et al.l ( 120071 ) at 
a center frequency of 1518 MHz. These pulsars were PSRs J1049— 5833, J1502— 5653, and 
J1525— 5417 (see their Table 1). We compared our NF values as a check for consistency. 
Two of the archival Parkes observations that we used for this comparison were also taken 
at 1518 MHz, while the third (PSR J1049-5833) had a center frequency of 1318 MHz. The 
uncertainties in our measured NFs were estimated using N~ 1 / 2 , where N is the number of 
su bintegrations (pulse s) used. This is a slightly more conservative estimate than that used 
by IWang et al.l (120071 ). where the number of null subintegrations n p and total number of 
subintegrations N were used in the e xpression ^/nZ/N . In all three test cases we obtained 
values that were consistent with the IWang et al.l (120071 ) values considering our respective 
uncertainties and that different data sets were used. This NF comparison is presented in 
Table [2j As mentioned below, the NF measured for PSR J0529— 6652 may be influenced by 
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the much larger distance to the pulsar than to the test pulsars, possibly making many of the 
pulses from PSR J0529— 6652 unseen. 



Results and Discussion 



3.1. Giant Pulses from PSR J0529-6652 and the Nulling Fraction 



Individual pulses from PSR J0529— 6652 were clearly detected in the 4195.5 s subset 
of the observation used in the analysis (Fig. [3]). Three of these pulses, shown in the lower 
part of Fig. [3J, were constructed using the dedispersed data. In these three plots, the 
intensity (shown on the vertical axis) was calculated by equating the rms of the dedispersed 
time series to the noise level as determined by the radiometer equation (see, e.g., Table 1 of 
Burke-Spolaor et al. 2011). The pulse flux density was calculated from the rms and the pulse 
signal-to- noise ratio (S/N), which ranged from 8 to 11 for these three pulses. The resulting 
flux densities ranged from S ~ 420 to 575 mJy at 1400 MHz. Using the definition of the 
pseudoluminosi ty (e.g., Lorim er Sz Kramer 2005), L = Sd 2 , and assuming a LMC distance 
of d ~ 50 kpc (|Koerwerll2009l ). this corresponds to 1400 MHz peak luminosities of between 
1050 and 1440 Jy kpc 2 . All three pulses have widths of ~ 20 ms, corresponding to about 40 
time samples, which is comparable to the width of the integrated pulse profile shown in Fig. 
[21 This is much larger than either the dispersion smearing within the frequency channels 
(~ 0.15 ms) or the expected pulse scattering time (~ 10~ 4 ms) from Galacti c plasma at 1400 
MHz, as estimated from the NE2001 electron model (jCordes &: Lazioll2002l ). This indicates 
that the pulses are not giant micropulses or very narrow single pulses like those seen for the 



Crab and Vela pulsars an d for PSR B1937+21 (jCognard et al. 



19961 : IJohnston et all 12001 



Johnston fc Romanil 120041 ) . The three pulses shown also occur at the same pulse phase as 



the peak of the integrated profile. 

A comparison of pulse intensities in the sample to the average pulse intensity indicates 
that some pulses have intensities > 20 times the average. This can be seen in Fig. HI which 
shows the pulse intensity histogram. These pulses may be consi dered giant pulses, which are 



generally defined as having amp litudes > 10 times the average (IJohnston et al.ll200ll ; iKnight 



20061 ; iKaruppusamy et al.ll201ll ). However, it is not clear from the small number of high- 



intensity pulses shown in the histogram in Fig. H] whether the giant pulses follow a power-law 
distribution in intensity, as expected for classical giant pulses (e.g., Argyle & Gower 1972), 
or possibly a lognormal distribution, as in the case of PSR B0656+14 (Weltevrede et al. 
2006b; see discussion below). A fit to the data is unable to distinguish between these two 
cases. Future observations should provide better statistics and should be able to resolve 
this question. There is also no evidence of any association with a high-energy emission 
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mechanism for PSR J0529 — 6652, which is thought to be a feature of classical giant pulses 
( [Johnston &: Romanil 12004 ) . We checked this with a high-resolution XMM-Newton map of 
the region (Filipovic, priv. comm.; see also Bozzetto et al. 2012) and by folding Fermi Large 
Area Telescope (LAT) data using the pulsar's ephemeris to look for gamma-ray pulsations. 
No emission was detected from the pulsar in either case. 



Another test for the presence of giant pul s es is a n R-parameter, defined by lJohnston et al 



(120011 ) and described by lBurke-Spolaor et al.l ( 120121 ). A value R is determined for each profile 
bin by taking the difference between the maximum value and the mean value of the samples 
for that bin, which is then divided by the rms of the values for that bin. A subsequent 
comparison of the R-value for each bin with the rms of the R-values for the off-pulse bins 
can be used as a t hreshold test for the presen ce of sparse modulated emission, such as giant 
pulses. Following iBurke-Spolaor et al.l ( 120121 ). we use a significance threshold defined as R 
minus the off-pulse mean of the R values, divided by the standard deviation of R in the 
off-pulse window. If this quantity is greater than 4 for any phase bin, then giant pulses are 
considered to be present. For three of the four on-pulse bins for PSR J0529— 6652, R was 
larger than 6 and nowhere else in the profile was it above 4. This supports our conclusion 
that giant pulses from PSR J0529— 6652 are present and that they do not significantly lead 
or trail the main pulse. 

As seen in Fig. [3J there are at least 15 distinct pulses th at are visible by eye that wer e 
also detected in the single pulse search detection algorithm ( jCordes fc McLaughlin! 120031 ). 
This corresponds to one detectable pulse every ~ 4.5 min on average, though it is evident 
that the pulses are not evenly spaced. For comparison, we conside r two young;, ca n onica l 
giant radio pulse emitters: the Crab pulsar and PSR B0540— 69. ICrawford et al.l (120051 ) 
used the same Parkes observing system as used here (but with a different sampling rate) 



to ob serve the X-ray pulsar PSR J0537-6910 in the LMC at 1400 MHz (Marshall et al. 



19981 ). They calculated that if the Crab pulsar were located in the LMC, one pulse would be 
expected to be detectable with this system every ~ 20 min. In the case of PSR B0540— 69, 
which is in the LMC, one giant pulse should be detectable every ~ 30 min@ Thus, the rate of 
detectable pulses per unit time from PSR J0529— 6652 with this system is almost an order of 
magnitude greater than the expected rates for either the Crab pulsar (if it were in the LMC) 
or PSR B0540— 69. Moreover, the spin period of PSR J0529— 6652 (976 ms) is more than an 
order of magnitude larger than either the Crab pulsar (33 ms) or PSR B0540— 69 (50 ms). 
Therefore, the likelihood that any one pulse from PSR J0529— 6652 would be detectable as 



6 This expected detection rate for PSR B0540— 69 is somewhat larger than the actual detection of only 
one pulse in a recent 1.4 hr test observation that had comparable sensitivity. This observation was taken as 
part of new multibeam pulsar survey of the LMC with Parkes. 
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a giant pulse with this system (or, alternatively, the fraction of emitted pulses detectable as 
single pulses with this system) is about two orders of magnitude greater than for either the 
Crab (if it were in the LMC) or PSR B0540— 69. Unlike these two pulsars, however, PSR 
J0529— 6652 is physically unremarkable, making this feature of its emission quite unusual 
and unexpected. 

The pulse c haracteristics from PSR J0529— 6652 appe ar to be more similar to those from 
PSR B0656+14 rtWeltevrede et al.ll2006ri iTao et alJlioii ). The pulses fro m PSR B0656+14 



follow a lognormal distribution rather than a power-law (e.g., Fig. 1 of IWeltevrede et al. 



( I2006bl ). which might also be the case for PSR J0529— 6652 (see Fig. HJ), they are not associ- 
ated with any high-energy emission, and they are not narrow micropulses. PSR B0656+14, 
like PSR J0529— 6652, also does not have a large magnetic field at the light-cyl inder radius, 
which is a feature of known classical giant pulse emitters ( ICognard et al.l 119961 ) . These fea- 
tures make both PSR J0529— 6652 and PSR B0656+14 different from classical giant pulse 
emitters. As discussed below, the large modulation indices measured for PSRs B0656+14 
and J0529— 6652 and the proximity of both pulsars to the bulk of the unrecycled pulsar 
population on the period/period derivative diagram also suggest a similarity between these 
two sources. 

We measured NFs for PSR J0529— 6652 and our three test pulsars. These values are 
presented in Table [TJ Both PSRs J0437— 4715 and J1359— 6038 have essentially no nulling 
near 1400 MHz (NF ~ 0). The NF for PSR J0529-6652 at 1390 MHz is 83.3 ± 1.5%, 
indicating that either radio emission is not present for the majority of the pulse periods or 
that most pulses are simply too weak to be distinguishable from noise. This latter point is 
an important possibility given the very large distance to PSR J0529— 6652. In all four cases 
the measured NFs appear qualitatively consistent with the pulse stacks shown in Fig. [TJ 



3.2. The Modulation Index of PSR J0529-6652 

For the modulation index calculation, we first estimated the contribution to the modula- 
tion from the Galactic ISM for each of the pulsars, and we corrected the measured modulatio n 



indices for this to obtain the value intrinsic to each pulsar. As outlined in lJenet &: Gill ( 120031 ). 
the measured (m), intrinsic (mj), and ISM-produced {ttiism) modulation indices are related 
according to: 



(m 2 + l) = K 2 + l)K SM + l). 
mjsM can be estimated using the relationship ( Uenet et al 



2001 



Jenet & Gil 



(4) 

2003h 
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miSM 



t)B/8v 



A/2- 



(5) 



where B is the observing bandwidth, bv is the characterist ic ISM diffractive sc intillation 



bandwidth, and 77 is a coefficient ranging from 0.1 to 0.2 (ICordes et al. 



199 01). 8y was 



determined at 1400 MHz for each pulsar using the NE2001 model of ICordes fc Laziol (120021 ) 
and the pulsar's catalog DM. For PSR J0529— 6652, we used the maximum Galactic DM 
contribution along its line of sight according to the NE2001 model (52 pc cm -3 ). The 
uncertainty i n itits m was calculated using a range of 77 of 0.1 to 0.2. Following the method of 
Jenet fc Gill ( 120031 ). the uncertainty in the intrinsic modulation index was determined by 
using the larger of the measurement uncertainties in m and ttiism- The results are presented 
in Table [TJ 

It should be noted that in the case of PSR J0529— 6652, uiism does not include any 
contribution from propagation through the LMC, where the pulsar resides. There is no model 
that we can use to estimate the LMC contribution to the modulation, but observations of 
other known pulsars in the MCs do not sho w any of the extrem e variability and null-lik e 
behavior that is seen for PSR J0529-6652 JCrawford et al.lboOll ; iManchester et all bood ) . 
We conclude from this that the contribution to the modulation from the LMC itself is likely 
to be minimal and does not significantly affect our results. 

After correcting for the estimated Galactic ISM contribution, we derived an intrinsic 
modulation index value for PSR J0529— 6652 for each profile phase bin. Fig. [5] shows the 
modulation indices for the on-pulse profile bins overlaid with the profile intensity. We selected 
the bin with the highest S/N (the profile peak) for the analysis , whic h also had the most 
precise modulation index measurement. As stated in I Jenet fc Gill ( 120031 ). this is the bin that 
is least likely to be affected by outlying conal emission components, which can introduce 
additional variability and can complicate the interpretation of the modulation index as a 
feature of the polar cap physics. This bin also had the lowest value of mj(0), as expected 
(see Jenet Sz Gil 2003), and from this we derived an intrinsic modulation index of 4.07 ±0.29 
for PSR J0529— 6652 (Table [TJ). The modulation indices measured for our three test pulsars 
using the same method are also presented in Table [TJ As expected, these values are all small 
(rrii < 1) and are consistent with the range seen for most radio pulsars (see, e.g., Fig. |6|). 

Fig. [TJ shows the pulse stack for PSR J0529— 6652, which clearly shows its variability. 
As stated above, the LMC itself is likely not a major contributor to this variability. The 
diffraction s cintillation time scale at 1400 MHz for PSR J0529— 6652 from Galactic plasma 
is ~ 1500 s (ICordes fc Laziol 120021 ). significantly longer than the minute time scales (or less) 
that seem to be present qualitatively in the pulse stack. The scintillation bandwidth is also 
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much smaller than the observing bandwidth (see Tabled]). This suggests that the variability 
of the emission is intrinsic to the pulsar and is consistent with the large intrinsic modulation 
index measured. 

Fig. [6] shows the intr insic modulation index vs. the complexity parameter derived using 



the 


Gil & Sendvk 


(2000) 


Weltevrede et al. 


'2006a) 



d (|2000 ) model for several samples of pulsars. 174 pulsars measured by 



that their Crab measurement was excluded i n their figure and is a lso no t included here). In 



addition, 102 of the 103 pulsars measured by lBurke-Spolaor et al.l (120121 ) for which there was 



also a cataloged P are plotted as small red dots, but no error bars are given for these values. 
Our three test pulsars are plotted as small black dots with error bars, and our measurement 
for PSR J0529— 6652 is indicated by the large diamond (see also Table [Q. All of these 
measurements were taken near 1400 MHz. Although PSR J0529— 6652 has a complexity 
parameter and physical characteristics that are typical of the rest of the sample (and the 
overall radio pulsar population in general; see Fig. [7]), it nevertheless has a modulation 
index that l i es we ll above this sample and the sample of 12 pulsars shown in Fig. 1 of 



Jenet fc Gill (120031 ) . Only two of the pulsars measured in any of these samples (apart from 



PSR J0529-6652) have m > 2. 

The large modulat ion index for PSR J052 9— 6652 is closer to the ranges and limits 
presented in Table 4 of IWeltevrede et al.l (120111 ) for four members of different neutron star 
classes (these are plotted in Fig. E]as large squares). These four sources are the Crab pulsar 
(see above), the radio magnetar XTE J1810-197, the RRAT J1819-1458, and the RRAT- 
like pulsar PSR B0656+14. All of the measurements of these sources were taken near 1400 
MHz, except for PSR B0656+14, which was taken at 327 MHz. Like PSR J0529-6652, 
they all have modulation indices that are significantly larger than the values that have been 
measured for typical radio pulsars. However, with the possible exception of PSR B0656+14, 
none lies near the center of the unrecycled radio pulsar population on the period/period 
derivative diagram (Fig. [7}. Thus, only PSR J0529-6652 (and possibly PSR B0656+14) 
are typical pulsars in this sense. In any case, the core-like characteristics of the emission from 
PSR J0529— 6652 (suggested by its simple and narrow profile morphology and its symmet- 
ric sign change in circular polarization), plus its unremarkable spin and physical parameters 
(Crawford et al. 2001; see al so Fig. 6), suggest th at this pulsar is more useful for a test of the 
correlation predicted by the iGil fc Sendykl ( 120001 ) emission model than these other neutron 
stars, all of which are physically unusual in different ways. We conclud e that the larg e mo d- 
ulation index measured for PSR J0529— 6652 is not consistent with the lGil fc Sendykl (120001 ) 
emission model, which predicts a small modulation index for a relat ively large compl exity 
parameter. We note also that the radio magnetar PSR J1622— 4950 ( Levin et al. 201 cl) has 
a relatively small measured modulation index (m ~ 0.5 measured by iBurke-Spolaor et al. 
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( 120121 ) at 1.4 GHz, and m ~ 1.7 measured by lLevin et al.l (120121 ) at 3.1 GHz). This indicates 
that not all magnetars exhibit extreme radio modulation like XTE J1810— 197, just as PSR 
J0529— 6652 illustrates that not all typical radio pulsars have low modulation. 



4. Conclusions 

We have studied the variability of PSR J0529— 6652, a luminous radio pulsar in the LMC, 
using observations taken at 1390 MHz with the Parkes 64-m telescope. The pulsar emits 
detectable single radio pulses that can be classified as giant pulses, making this the second 
known pulsar in the LMC (after PSR B0540— 69) to do so. These pulses are comparable in 
width to and occur at the same phase as the integrated profile peak, suggesting that they are 
not giant micropulses. The characteristics of the giant pulses appear to be closer to those seen 
from PSR B0656+14 than from classical giant pulse emitters, such as the Crab pulsar. The 
fraction of the pulses emitted by PSR J0529— 6652 that are individually detectable as single 
pulses at this frequency is two orders of magnitude greater than either PSR B0540— 69 or the 
Crab pulsar (if the latter were located in the LMC). PSR J0529— 6652 also appears to exhibit 
nulling behavior, and we have measured a NF of 83.3 ± 1.5% for the pulsar. Pulsed radio 
emission is either not present for the majority of pulse periods, or it is too weak in most cases 
to be individually detectable above the noise. Given the large distance to the pulsar, this 
is a possibility that must be considered. Our measured intrinsic modulation index for PSR 
J0529— 6652 is 4.07 ± 0.29, which is significantly larger than the values previously measured 
for typical radio pulsars w ith similar spin and phys ical characteristics. It is comparable to the 



larger values presented by lWeltevrede et al.l (120111 ) for several members of other neutron star 



source classes that are known to be radio variable. Unlike these other sources, however, PSR 
J0529— 6652 has spin and physical characteristics that are typical of the un recycled radio 



pulsar population. These features make this pulsar useful as a test of the iGil fc Sendyk 



( 120001 ) emission model, and the large modulation index measured for PSR J0529— 6652 does 
not appear to be consistent with the model prediction. This conclusion depends to some 
degree on the assumption that PSR J0529— 6652 is exhibiting core emission, as seems to be 
the case given its simple profile morphology, narrow profile width, and previously measured 
polarization characteristics at 600 MHz. 



We thank M. Filipovic for providing us with an X-ray image of the region containing 
PSR J0529— 6652. The Parkes radio telescope is part of the Australia Telescope, which is 
funded by the Commonwealth of Australia for operation as a National Facility managed 
by CSIRO. This work was supported in part by grants from Research Corporation, the 
Mount Cuba Astronomical Foundation, the National Radio Astronomy Observatory, and 
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Fig. 1.— Pulse stacks for PSR J0529-6652 and three bright test pulsars (PSRs J0437-4715, 
J0536— 7543, and J1359— 6038). Each pulse stack has 128 pulse phase bins (horizontal axis) 
but a different number of consecutive pulses (vertical axis). Table [1] lists the observing 
parameters and ATNF catalog properties for the four pulsars. The contrast in each plot has 
been adjusted to best illuminate any variability in the pulses. All observations were taken 
near 1400 MHz, and the data are largely free of RFI. The variability is clearly evident in the 
PSR J0529— 6652 pulse stack, consistent with its large measured NF (Tabled]). 
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Fig. 2. — Normalized mean intensity profile for PSR J0529— 6652 from the addition of 4299 
consecutive pulses from a 1390 MHz Parkes observation. There are 128 phase bins in the 
profile, and the profile has unity peak value and an off-pulse mean of zero. The mean pulse 
profile is narrow and uncomplicated, with a width of 3-4 bins (~ 25 ms, or ~ 3% of the 
pulse period), and has no obvious additional or outlying components. These features, plus 
the polarization characteristics measured at 600 MHz by I Cost a et al.1 (Il99ll ). suggest that 
PSR J0529— 6652 is likely to be exhibiting core emission. 
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Fig. 3.— Single pulse detections of PSR J0529-6652 at 1390 MHz from 4195.5 s of inte- 
gration time, or 4299 pulses. The data here are the same data that were used to make the 
pulse stack and to measure the modulation index and NF (see Table d]). The top frame 
shows pulse strength as a function of both DM and time, with pulse events with S/N > 5.5 
shown. As expected, the pulses occur most strongly near the pulsar's DM of ~ 100 pc cm~ 3 . 
The lower frames show three detectable pulses from this observation plotted as intensity vs. 
time after dedispersion was applied. In all three cases, the pulses have a width of ~ 20 ms 
(~ 2% of the pulse period), or roughly 40 samples. This is comparable to the width of the 
integrated pulse profile (see Figs. [2] and [5]). Since the pulses are not dispersion or scatter 
broadened, this suggests that we are seeing the intrinsic widths and that the pulses are not 
giant micropulses. These three pulses also occur at the same pulse phase as the integrated 
profile. PSR J0529-6652 is the second pulsar in the LMC (after PSR B0540-69) known 
to emit detectable single radio pulses. The pulses shown here illustrate the high degree of 
amplitude variability for the pulsar, which is confirmed by the large measured modulation 
index and NF. 
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Fig. 4. — Histogram of on-pulse intensities for PSR J0529— 6652 from the subset of 4299 
pulses used in the analysis (solid line). The corresponding off-pulse intensities calculated 
using the same number of off-pulse phase bins for each pulse are also shown (dashed line). 
Both histograms have been normalized to the mean on-pulse intensity value. These his- 
tograms were used in the calculation of the NF for PSR J0529— 6652. There is an excess 
of pulses with large amplitudes, extending well beyond the noise limit, indicating that PSR 
J0529— 6652 is a giant pulse emitter. We are unable to distinguish between a power-law and 
lognormal distribution for the giant pulses owing to the relatively small number of pulses 
detected. 
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Fig. 5. — Phase-resolved intrinsic modulation index values measured for PSR J0529— 6652 
(squares) overlaid with corresponding mean intensity values (dash-dotted line). The intensity 
values have been arbitrarily scaled for display purposes. Only the pulse phase bins near the 
pulse peak are shown, and only the on-pulse bins have enough signal for reliable modulation 
index measurements. These modulation index values have already been corrected for the 
estimated contribution from fluctuations from the Galactic ISM. The minimum and most 
precise value of mi = 4.07 ± 0.29 is seen at the profile peak, and this is the value used in our 
analysis (see, e.g., Jenet & Gil 2003). 
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Fig. 6. — In t rinsic modulation index vs. complexity parameter determined from the 
Gil fc Sendvkl feoooh model is shown for several samples of pulsars. All of the 174 pul- 



sars measured by IWeltevrede et al.l ( )2006al ) (and presented in their Table 2 and Fig. 10) are 
presented as small blue circles with error bars. Also plotted are the three measurements of 
o ur test pulsars (small black circles with errors). The 102 (out of 103) pulsars from Table 1 
of iBurke-Spolaor et al.l ( 120121 ) with cataloged P values are plotted with red circles (no errors 
were reported). Ou r measurement for PSR J0529— 6652 (large diamond) and the measure- 
ments presented by IWeltevrede et al.l (120111 ) in their Table 4 for four members of different, 
more variable neutron star classes are also shown (large squares). All data shown here were 
taken near 1400 MHz, except for PSR B0656+14, which was at 327 MHz. PSR J0529-6652 
has physical parameters that are typical of the unrecycled radio pulsar population (see Fig. 
[7]), unlike the other labeled neutron stars (except possibly PSR B0656+14). 
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Fig. 7. — Period der ivative vs. period for radio pulsars from the ATNF catalog 



([Manchester et al.ll2005[ ) (small dots). Not all pulsars are shown (e.g., recycled millisecond 
pulsars are beyond the plot limits). Also shown are the four neut r on sta rs with modula- 
tion index measurements presented in Table 4 of IWeltevrede et al.l (120 111 ) (squares). PSR 
J0529— 6652 is shown as the large diamond. Compared to the other four labeled neutron 
stars, PSR J0529— 6652 lies closer to the center of the unrecycled radio pulsar population 
and has spin characteristics that are more typical of this population. 
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Table 1. Properties, Observing Parameters, and Measured Quantities for Four Pulsars 



Pulsar 


J0529-6652 


J0437-4715 


J0536-7543 


J1359-6038 


Topocentric period, P (s) a 


0.97573571876 


0.00575786 


1.24585559629 


0.12751295 


Dispersion measure, DM (pc cm~ 3 ) a 


103.2 


2.6 


17.5 


293.7 


Number of consecutive pulses used, iV 


4299 


20839 


240 


1943 


Number of profile bins 


128 


128 


128 


128 


Center observing frequency (MHz) 


1390 


1390 


1390 


1374 


Sampling time (/is) 


500 


80 


250 


500 


Observing bandwidth, B (MHz) 


256 


256 


256 


288 


Scintillation bandwidth, 8v (MHz) b 


1.51 


318 


22 


~ 


Modulation index for the ISM, tuism 


0.20 ±0.03 


0.95 ±0.02 


0.61 ±0.07 


~ 


Measured modulation index, m d 


4.15 ±0.29 


1.00 ±0.02 


1.40 ±0.22 


0.33 ±0.03 


Intrinsic modulation index, m, c 


4.07 ±0.29 


0.23 ±0.02 


1.08 ±0.22 


0.33 ±0.03 


Measured nulling fraction (NF) f 


83.3 ± 1.5% 


0.0 ±0.7% 


32.5 ±6.5% 


0.1 ±2.3% 



Note. — All observations were taken with the center beam of the multibeam receiver and 512 x 0.5 MHz 
channel analog filter bank system at the Par kes telescope, except for PSR J1359— 6038, where 96 x 3 MHz 
channels were used ([Manchester et al.ll2001l ). 

a Obtained or derived from the ATNF pulsar catalog ([Manchester et alJl2005l ). 



b Estimated from the NE2001 Galactic electron model of ICordes fc Laziol ([2002). 

c Estimated contribution to the modulation from propagation through the Galactic ISM. 

d Measur ed from the norma lized pulse stack. The lowest value of m((f>) among the on-pulse bins was 
chosen (see lJenet fc Gil d2003h and the text for justification). 

c Obtained after correction for the estimated Galactic ISM contribution. 

f The uncertainty in NF was determined by iV -1 / 2 , where N is the number of pulses used. 
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Table 2. Nulling Fractions for Three Known Nulling Pulsars 



Pulsar 


Wans et al. (20Q7) a 


This Work 


iV b 


J1049-5833 


47 ± 3% 


33 ± 35% 


8 


J1502-5653 


93 ± 4% 


70 ± 9% 


120 


J1525-5417 


16 ± 5% 


26 ± 5% 


361 



Note. — For all three pulsars, different sets of Parkes 
data were u s ed for our measurements than were used by 



Wang et al.1 (|2007l ). All of the data used in both sets of 



measurements had a center frequency of 1518 MHz, except 
for our measurement of PSR J1049— 5833, which had a 
center frequency of 1318 MHz. 

a Values from Table 1 of 1 Wang et all |2007h . 

b Number of subintegrations used in our NF measure- 
ments. The uncertainty in the NF was determined by 

JV-l/2. 



